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Intensity  Scintillations  of  an  EM  Wave  in 
Extremely  Strong  Turbulence 


I 

I: 


j 


1.  INTRODUCTION 

1-4 

Recent  theoretical  research  on  the  propagation  of  electromagnetic  waves 
in  extended  strong  turbulence  indicates  that  the  electric  field  obeys  gaussian  sta- 
tistics in  the  limit  when  the  strength  of  the  turbulence  becomes  infinite.  In  fact, 

5 

it  has  been  shown  that  theoretical  solutions,  which  have  as  their  basis  the  assump- 
tion of  gaussian  statistics  for  extremely  strong  turbulence,  give  excellent  agree- 
ment with  measured  data  for  the  intensity  fluctuations  in  the  case  of  moderately 
strong  turbulence.  Therefore,  even  though  measured  data  indicate  that  the  field 
statistics  are  nearly  log-normal  in  moderately  strong  turbulence,  there  is  good 
reason  to  believe  that  the  field  statistics  are  gaussian  when  the  turbulence  is 
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extremely  strong.  Unfortunately,  there  are,  at  present,  no  measurements  of  the 
field  statistics  in  extremely  strong  turbulence. 

In  this  paper,  we  shall  assume  that  the  field  statistics  of  a point  source 
approach  gaussian  in  extremely  strong  turbulence.  By  using  that  assumption,  we 
shall  then  calculate  expressions  for  the  averaged  intensity,  intensity  correlation 
and  scintillation  index  for  a source  of  arbitrary  shape;  the  source  may  be  either 
coherent  or  incoherent.  We  find  that  the  scintillation  index  always  approaches 
unity.  This  result  appears  to  disagree  with  some  recent  calculations  by  Rumsey6 7 
who  found  that  the  scintillation  index  for  an  Incoherent  source  approaches  zero  in 
extremely  strong  turbulence.  The  reason  for  this  apparent  discrepancy  is  inves- 
tigated. 


2.  COHERENT  SOURCE 


Let  us  first  consider  the  scintillations  in  the  received  signal  from  a coherent 

source,  such  as  a laser.  By  using  the  extended  Huygens -Fresnel  principle,  we 

can  show  that  the  field  distribution  e(p)  in  the  plane  x is  related  to  the  source  dis- 
7 — 

tribution  e^jj^)  via 


e0(£l)  exp 


k 2 

^ +^£*£1> 


(1) 


where  k is  the  signal  wavenumber  and  is  the  additional  complex  phase  due 

to  turbulence  of  a spherical  wave  propagating  from  (0,j>jj)  to  (x,£>.  We  now  desire 
to  obtain  the  long-term  averaged  intensity  (1(g))  and  the  long-term  averaged  rnten- 
sity  correlation  (%)I(£')).  By  long-term  averaged,  we  mean  that  the  average  is 
taken  over  times  long  compared  with  the  coherence  time  of  the  turbulence.  For 
example,  for  a beam  of  diameter  D we  take  the  average  over  times  long  compared 
with  D/V,  where  V is  the  transverse  velocity  of  the  turbulent  eddies.  By  using 
Eq.  (1)  we  find  i'  straightforward  to  show  that 


6.  Rumsey,  V.  (1975)  Scintillations  due  to  an  extended  source  in  extended  strong 

turbulence,  Digest  of  the  1975  Annual  USNC/URSI  Meeting.  Nat'l  Academy 
of  Sciences,  Wash.,  D. C.  20418,  p 246. 

7.  Feizulin,  Z. , and  Kravtsov,  Y.  (1967)  Broadening  of  a laser  beam  in  a turbu- 

lent medium,  Radlophys.  Quantum  Electron  10:33-35. 
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(%))  a (e(g_)  e*(£)) 

= (2^)  / •**  /d2Pld2p2eo(£l)eo(£2) 

“oo  “oo 

• exp  {*££  ‘ (pi  -£2)  + i(Pi  "P^J 

. (exp[0(£,£ x)  4 i//*  <£,£2)]> 


and 


,4  co 


<I(£)I(£,))  = J...Jd2Pl  d2p2  d2p3  d2p4  eo(£l)  e*(£2)  eo(£3)  e*(£4) 

exp  j 1 x £ ' fel  “£2^  + 1 x£'  ‘ ^£3  “£4^  +-2x  *P1  ~ p2  + p3  ’ 


• (expf^p.p^  + ^*(p,£2>  + <M£',£3)  + ^*(£',£4)]  ) 


(3) 


Let  us  now  employ  the  ergodic  hypothesis  so  that  the  time  average  is  replaced 
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by  an  ensemble  average.  In  this  case  it  can  be  shown  that 

(exp[0(£f£l)  + tf'*(£',£2)])  = exP  |“  \ di(£  £1  "£2)|  (4) 

O 

where  for  a Kolmogorov  spectrum  for  the  index  of  refraction  fluctuations,  we  have 

D1(|,r1)  = k 2f  dx1  C2(x')  g(/3)  (5) 

*0 


8.  Kon,  A.,  and  Feizulin,  Z.  (1970)  Fluctuations  in  the  parameters  of  spherical 
waves  propagating  in  a turbulent  atmosphere,  Radiopftys.  Quantum  Electron. 
13:51-53.  
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with 


¥ 


i 


„ _ _ X'  /,  X'  ' 

0=  It  n 

) ■ 

(6) 

g(/3)  ~ 1.  8 a2  kJ/3 

if 

0 « 1/Km  . 

(7a) 

g(0)  ~ 2.92  /35/3 

if 

(7b) 
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and  C_  is  the  index  of  refraction  structure  constant  and  K,  - 5.91/f  , 
n m ' 0 

the  inner  scale  size  of  the  turbulence. 

where  i is 
0 

We  now  assume  that  the  turbulence  is  sufficiently  strong  so  that  the  function 
■ exp  is  a gaussian  random  variable.  As  discussed  earlier,  the 

general  theoretical  solutions  of  Prokhorov  et  al3  and  Fante,  * as  well  as  the  phys- 
ical model  of  de  Wolf,  1 suggest  that  this  should  be  the  case  for  sufficiently  large 
values  of  k7 ^ C3  xA*/*\  For  gaussian  point -source  statistics  we  have 

(exp x)  + £2)  + 0(£',£3)  + <£'.£4)]) 

= <exp[0(£,£!)  +/(£.£2)]><exP^(£'»£3)  +V'>K(£I,£4)]> 

+ (exp[0(g>i£1)  + 0*(£',£4))><exp(0(e',£3)  + 0*fe,£2)J> 

= exp  £-  i D1(0,£l  -£2)  - i -£4)J 

+ D1(£  -£',£l  -£4)  -I  -£,£3  -£g)j  . (8) 

If  we  now  use  Eqs.  (4)  and  (8)  in  Eqs.  (2)  and  (3)  we  get 


<!(£))  = 


d <>2  V-Sl*  eo— 2}  exP 


• exp 


k 

x 


P ■ (p  1 


(9) 
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and 


^ CO 

<I(£)I(£*)>  -■  (2^)  J“'Jd2P  1 d2p2  d2P4  eo(Pl}  e0(f2}  eo(£3)  *lbA) 


• exP  jljf.fi  • tel  -J>2)  + i|£'  ' (^3  '-£4)+fe(p5  -P2+P3  'p2) 

• j exp  ^|d1(0,£i  -£2)  -£4)J 

+ exp  -£',  £X  -£4)  -jDjfe'  -£,£3  -£2)jJ  . (10) 

Now  let  us  consider  the  scintillat  index 


m2 , <'2<p»_-  w.2.  , 

<i<p» 

Upon  setting  p'  - p and  substituting  Eqs.  (10)  and  (9)  into  Eq.  (11),  we  find  that 

m2  = 1 (12) 

2 

Therefore,  for  an  arbitrary  coherent  source  m — 1 as  the  streng  turbulence 

approaches  infinity.  This  result  is  a generalization  of  earlier  results  which 
4 2 

showed  that  m — 1 for  a plane  wave. 


3.  INCOHERENT  SOURCE 

We  now  consider  the  case  when  the  source  is  incoherent.  We  shall  assume  that 
the  coherence  time  Tg  of  the  source  is  much  smaller  than  the  coherence  time  T&  of 
the  atmospheric  turbulence.  Therefore,  we  shall  first  obtain  the  average  of  %) 
and  I(£)  • I(£')  over  times  which  are  long  compared  with  rg  but  short  compared 
with  For  an  incoherent  source  containing  very  many  coherence  areas  (the 

a 

coherence  area  is  of  linear  dimension  cr.  where  c is  the  c leed  of  light),  it  can  be 
9 s 

shown  that 


9.  Beran,  M.,  and  Parrent,  G.  (1964)  Theory  of  Partial  Coherence.  Englewood 
Cliffs,  N.J.,  Prentice  Hall,  Englewood  Cliffs,  N.J. 
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(13) 


5 k / eo^l}  eo(£2)  dt  - h^lU{£-l  '■£2) 


where 


Ts  <<  T <<:  Ta  1 


Also 


eofel>  eo(-£2}  eo(^3}  eo(p4)  = V-£l}  6(£l  ~&2}  6(£3  "£4) 

+ w h^3]  6(-fil  ■£4)  6{&2  "£3}  (14) 

where  Is  is  the  intensity  distribution  of  the  source.  Upon  using  Eqs,  (13)  and  (14) 
in  Eqs.  (2)  and  (3),  one  finds  it  easy  to  show  that  the  average  over  times  T such 
that  T « T « t„  is 

S a 


/ \2  00 

I(£>  = (ife)  ff  d2pi  W exp  + 


and 


I(£)I(e')  = I (p)  %')  + 


&)f[ 


d Pi  w 


exp 


i - (p  "£')  ■ Pj  + ^(p,^)  + ^ ’*'(£', £j) 


(15) 


(16) 


Equation  (16)  is  the  generalized  version  of  the  Brown-Twiss  formula  and  for  no 
turbulence  (i 1/  = 0)  it  reduces  to  the  usual  Brown-Twiss  result. 

Next  let  us  obtain  the  average  over  times  which  are  long  compared  with  T • 

If  we  again  assume  ergodicity,  so  that  we  replace  the  time  average  by  an  ensemble 
average,  and  then  use  Eqs.  (4)  and  (8)  to  calculate  (expfiMg,^)  + </'’!t(£,£1)J)  and 
(exptyfg,^)  + ^(g'.jjj)  + ^*(£,£2)  + i/'(£i,£2^))  we  obtain*  for  the  long-time  aver- 
ages of  1(g)  and  1(g)  I(g * ) 


jjc 

In  writing  Eqs.  (17)  and  (18)  we  have  assumed  that  the  measuring  apparatus  can 
essentially  measure  the  instantaneous  values  of  I;  that  is,  we  assume  that  the 
apparatus  response  time  is  short  compared  with  the  source  coherence  time. 
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<%»  (ife)  fjA\W  (l7) 

*oo 

W %'»  ■ (ife)4  / " / A 1 Aa  Vi'1  >.W  [l  * exp  <-  D1<£  •*’’*»  'l2”] 

* (s;)"  • JVp , A J is(£i>  i,fea>  «p  [l  I fe  -£'>  • fei 

-oo 

• j^exp  {-Dj(£  -£',0)}  + expJ-D^O^j  -£2))]  • (18) 

Equations  (17)  and  (18)  represent  the  averages  over  times  longer  than  v assum- 
ing that  the  field  statistics  for  a point  source  become  gaussian  (rather  than  og- 
nox  mal,  as  is  the  case  in  weak  turbulence);  in  the  limit  of  very  strong  turbulence 
it  is  appropriate  to  use  the  form  in  Eq.  (7a)  for  Dj  so  that 

v . | 2 (19) 

Djto.pj  - p2)  2 Ql  Pi  " P2'  * 


where 


1.8  KlJ3  k2 


j dx»  x'2  C2i 


(x1)  . 


distribution  of  the  form  Ig(pj)  * exp  (-p2/2a2)  it  is  readily  found 


If  we  use  a source 


from  Eqs.  (17)  and  (18)  that  (for  Qa  » 1) 

<>>  sM(lfe)  ■ 

<l2>=a,<1>2  jl*nW,  ■ 

By  a more  complex  calculation  we  also  find  that 
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and  so  on.  From  Eqs.  (20a)  to  (20c)  it  is  clear  that  the  statistics  of  the  beam  are 
not  necessarily  gaussian  even  though  we  have  assumed  that  the  point  source  sta- 
tistics are  purely  gaussian.  This  is  "lear  because  for  pure  gaussian  field  statis- 
tics we  would  have  (I2)  = 2!  (I)2,  (I3)  = 3!  (I)3,  and  so  on,  whereas  from  Eqs.  (20a) 
to  (20c)  we  see  that  this  is  the  case  only  in  the  limit  of  a plane  wave  (a  — co). 
For  the  case  of  finite  (but  large)  Q,  the  field  statistics  of  a finite  beam  are  almost, 
but  not  quite,  gaussian.  * 

o 

We  can  now  use  Eqs.  (20a)  and  (20b)  to  calculate  the  scintillation  index  m . 

We  find  that 

m2  = l + 2 — _ (21) 

1 + 4Qa^ 

o 

from  which  we  see  that  m - 1 as  Q - » . Note  that  Eq.  (21)  is  valid  only  for 
o 

Qa  » 1 because  of  our  assumption  of  very  strong  turbulence.  It  is  interesting  to 
note  that  the  results  in  Eqs.  (13)  to  (21)  can  also,  with  suitable  modifications,  be 
used  to  describe  the  scatter  of  a beam  from  a diffuse  surface  located  in  a strongly 
turbulent  medium. 

Before  concluding  this  section,  we  find  it  interesting  to  see  what  would  be 
obtained  if  we  could  not  measure  the  instantaneous  value  of  I( g_)  but  could  measure 
only  the  value  of  I(jd)  averaged  over  a time  greater  than  the  coherence  time  Tg  of 
the  star  (but  short  compared  with  ra).  As  before,  let  us  denote  this  value  by  I(£). 

If  we  took  the  long-time  average  of  I(£)  I(p 1 ) (as  opposed  to  the  average  of  the 
product  KjjIKjj')  of  the  instantaneous  intensities)  we  get  from  Eq.  (15) 

/ \4  00 

<%)%•)>  • [m)  f—f  ^i  ^2  W W i1  + exP  (-Di(p  -p,.pi~p2}]  ► 

(22) 

If  the  measurements  are  made  in  the  manner  described  here,  the  equivalent  scln- 
o —2  — 2 

filiation  index  m = (I)/(I)  - 1 is,  for  the  form  of  assumed  in  Eq.  (20) 

m2  = i — , (23) 

1 + 4Qa^ 

which  approaches  zero  as  Q — co.  A result  equivalent  to  Eq.  (23)  has  been  obtained 
by  Rumsey. 


*That  is,  for  gaussian  field  statistics  the  probability  distribution  of  the  intensity  is 
p(I)  =-pyexp  (-  which  leads  to  moments  (I^1)  = Nl  (I)1^. 
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Therefore,  the  result  we  obtain  for  the  scintillation  index  of  an  incoherent 

source  defends  on  how  the  intensity  measurements  are  made.  If  the  measurement 

apparatus  can  obtain  essentially  instantaneous  values  (that  is,  time  response  short 

2 

compared  with  source  coherence  time  rs)  of  I and  these  are  used  to  obtain  (1  ) and 
(I)  then  the  measured  scintillation  index  approaches  unity  in  extremely  strong 
turbulence.  However,  if  the  measurement  apparatus  can  measure  only  the  value  of 
I averaged  over  a time  long  compared  with  the  coherence  time  of  the  source,  we 
find  that  the  measured  scintillation  index  will  approach  zero. 
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